Elastic interaction between colloidal particles in confined nematic liquid crystals 
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The theory of elastic interaction of micron size axially symmetric colloidal particles immersed 
into confined nematic liquid crystal has been proposed. General formulas are obtained for the self 
energy of one colloidal particle and interaction energy between two particles in arbitrary confined 
NLC with strong anchoring condition on the bounding surface. Particular cases of dipole-dipole 
interaction in the homeotropic and planar nematic cell with thickness L are considered and found to 
be exponentially screened on far distances with decay length \dd = ^ • It is predicted that bounding 
surfaces in the planar cell crucially change the attraction and repulsion zones of usual dipole-dipole 
interaction. As well it is predicted that the decay length in quadrupolar interaction is two times 
smaller than for the dipolar case. 
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Colloidal particles in nematic liquid crystals (NLC) 
have attracted a great research interest during the last 
years. Anisotropic properties of the host fluid - liquid 
crystal give rise to a new class of colloidal anisotropic 
interactions that never occurs in isotropic hosts. The 
anisotropic interactions result in different structures of 
colloidal particles such as linear chains in inverted ne- 
matic emulsions (T], [1] , 2D crystals H and 2D hexagonal 
structures at nematic-air interface [3, Q . 

Study of anisotropic colloidal interactions has been 
made both experimentally and theoretically [l0|- 

16|. The first theoretical approach was developed in 
l|, llOj l with help of ansatz functions for the director and 
using multiple expansion in the far field area. Another 
approach [ill , [l^ gave possibility to find approximate 
solutions in terms of geometrical shape of particles. Re- 
cently authors of [T5l.ll6l| proposed a method for finding 
elastic interaction between colloids based on the fixing of 
director field on the surface of virtual sphere surrounding 
the real particle. The predicted dipole-dipole forces are 
three times weaker and quadrupole-quadrupole five times 
weaker than results of ^M,- On the other hand authors 
of recently have measured experimentally both inter- 
actions and found that experimental results are in accor- 
dance with Lubensky et al. prediction [l^j with about 
10% accuracy. This allows to justify assumptions of [l3| 
for spherical particles for intinite nematic liquid crystal. 
In this paper we suggest to generalize that approach for 
the case of the confined nematic liquid crystals as prac- 
tically always NLC has to be confined with walls, cells 
or containers. In a broader context, understanding the 
elasticity-mediated colloidal interactions in confined me- 
dia is of great importance not only in the field of regular 
thermotropic liquid crystals, but also for understanding 
interactions in more complex media with orientational or- 
der, for example, in solutions of DNA, f-actin and other 



biologically relevant molecules. Up to now almost all 
experimental studies did not take into account quanti- 
tatively confinement effects besides the article of Vilfan 
et al. [4]. In that paper authors have found exponen- 
tial screening effects for quadrupole-quadrupole interac- 
tion between spherical particles in homeotropic NLC cell. 
From our viewpoint there was only one theoretical ap- 
proach for description of colloidal particles in confined 
NLC performed in papers ^3,, il4,] . 

In this paper we propose the new approach for quan- 
titative description of the axial colloidal particles in con- 
fined NLC. This method enables to find self energy of 
one colloidal particle and interaction energy between 
two particles in arbitrary confined NLC with strong an- 
choring condition n^(s) = on the bounding surface. 
We apply general formulas to the particular cases of 
dipole-dipole and quadrupole-quadrupole interaction in 
the homeotropic cell and to the dipole-dipole interaction 
in the planar cell with thickness L. 

Consider axially symmetric particle of the size O.ly^m-^ 
IQjjLm which may carry topological defects such as hy- 
perbolic hedgehog, diclination ring or boojums. Director 
field far from the particle in the infinite LC has the form 
nx{r) ^ p^ + ic^,ny(r) = p^ + 3c^ with p and c be- 
ing dipole and quadrupole moment (we use another nota- 
tion for c with respect to the c in [l^], so that our c = |c ). 
It was found in [lo| that p = aa^,c = —f3a^ with a being 
the particle radius, and for instance a = 2.04,/? = 0.72 
for hyperbolic hedgehog configuration. In order to find 
energy of the system: particle(s) -|- LC it is necessary 
to introduce some effective functional Feff so that it's 
Euler-Lagran ge e quations should have the above solu- 
tions. In the [10] it was found that in the one constant 
approximation with Frank constant K the effective func- 
tional has the form: 
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Feff = K 



(1) 



which brings Euler-Lagrange equations: 



c'.a^C(x)] 



(2) 



where P(x) and C(x) are dipole- and quadrupole 
moment densities. For the infinite space 

the solution has the known form: ^^'(x) = 
Id^^'j^^[-dlPix') + d'^d'Mx')]. If we con- 
sider P(x) = pS{x.) and C(x) = cS{x) this really brings 
n.(r) =P^ + 3c^, ny{r) = + 3c||. 

In the case of confined nematic with the boundary con- 
ditions n^(s) = on the surface S the solution of EL 
equation has the form: 



n^(x) 



/ d3x'G(x, x') [-9;P(x') + d'^d'M^') 
Jv 



(3) 



where G is the Green function AxG(x, x') = — 47r(5(x— x') 
for X, x' e V and G(x, s) = for any s of the bound- 
ing surfaces. Consider N particles in the confined NLC, 
so P(x) = J2iPi^{^ - Xi) and G(x) = J^i^iSi^ - x^). 
Then substition © into F^ff brings: Feff = U^'^'^-f + 

U^r^teracUon ^j^g^.^ jjself ^ jjself ^ ^^^^ jjself 

interaction of the i-th particle with the bounding surfaces 



self 



u. 



self 
dd 



U^Q^ + Uqq . In general case the inter- 
action of the particle with bounding surfaces (self-energy 
part) takes the form: 

Ulf = ^2i,Kp^d,d'H{x,X)\^.=^' 



j-self 



where G(x,x') = + H(x,x') and Ax-ff(x,x') = 

|X X I 

(we excluded divergent part of self energy from j^jz^)- 
Interaction energy [/»»*eraction ^ j^^^- U^f. Here U^f 
is the interaction energy between i and j particles: 

Udd + UdQ + Uqq: 



J Tint 



Udd - -^TTKpp'dM{x^,x:A 



(5) 



UdQ = ~AnK {pc'a^9;9;G(x„x;) +p'ca;9^a,G(x„x;)} 



Uqq = -4^W9,aia^a;G(x„x;.) 

Here unprimed quantities are used for particle i and 
primed for particle j. Formulas ^ and ^ represent 
general expressions for the self energy of one particle 
(energy of interaction with the walls) and interparticle 
elastic interactions in the arbitrary confined NLC with 
strong anchoring conditions n^^is) = on the bounding 
surfaces. Below we will apply these expressions for par- 
ticular cases of the nematic cell with homeotropic and 
planar configurations. 



A. Interaction in the homeotropic cell with width L 

Green function in this case has the form : 



u 
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dQ 



47rKpcd^d'^d',H{x,,x'^) 



(4) 
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Here heights z,z', horizontal projections p<,p> and 
Im, Km are modified Bessel functions. Then using of (O 
brings dipole-dipole interaction in the cell : 



U. 



dd,hora 



IGnKpp' 
Z3 



> [mr) sm — — 



WKZ riTTZ ,nTTp 
sin——Ko[——) 



(7) 

with p being the horizontal projection of the distance 
between particles. Similar quadrupole-quadrupole inter- 



action takes the form: 



U, 



16ttKcc' -t-^ 4 mrz mrz' ^mrp 
2_^{n'K) cos—^cos—^Kq{—^) 

n=l 

(8) 

When both particles are located in the center 
of the cell z = z' = f we have U^^^^^^^ = 

^Tl^Er=i,o..("^)'^o(^) (see black thick line on 
the Fig. [T]). In the hmit of small distance p <^ L be- 
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FIG. 1. (Color online) Log-log plots of the interaction po- 
tential in kT units as a function of the rescaled interparticle 
distance p/ L. Here particle's radius a — 2.2/im,cell thick- 
ness L = 7iim,K = 7pN,p = p' = 2.04a^c = c' = 0.2a^. 
Blue thick line is quadrupole potential in homeotropic cell (it 
approximates experimental data from very well), dashed 
thick line is it's power-law asymptotics oc All thin lines 

are dipole potentials in planar cell. Coming anticlockwise 
first purple line is attraction along the direction ip = 40°, 
second brown line is along = 20°, third red line is along 
z axis (/p = and the last green thin line is repulsion along 
if — ^. Black thick line is dipole repulsion in homeotropic 
cell from Thin dashed line is the power-law asymptotics 

Tj 4nKp^ 



tween particles it has asymptotic 



^ dd.honi "f^ 

that is in agreement with standard formula for the usual 
dipole-dipole interaction Udd = ^^^^(1 - icos'^e) for 
9 = ^. From the Fig. [1] it is clearly seen that power- 



law behaviour U 



is valid to the p — 1.2L. For 



larger distances p > 1.2L exponential decay takes place 
Udd.r^oo = IGnKpp' 1/2Lp i^ave decay length for 

dipole-dipole interaction Add = ^■ 

When both particles are located in the center of the cell 
z — z' — ^ we have quadrupole interaction C/gg /jo„j = 

^T^Er=2.e.en("^)'^o(^)(see thick bluc line on 
the Fig. [T]) . This coinsides with the result of [13] if we 
take r there to be equal T = 2ttKc — —2f3-KKa^ . Let's 
emphasize that in [14| the F remains unknown quantity. 
In the limit of small distance p <^ L between particles 
it has asymptotics U[ 



ment with standard formula for the usual quadrupole- 
quadrupole interaction Uqq = AttKcc' ^-^OcosH+wbcos^e 
for 6* = ^. This power- lar behaviour is valid to the 



QQJiom 



distance p = 0.8L. For larger distance p > 0.8L 
crossover to the exponential decay occurs C/QQ,r=>oo — 

8t:Kcc'{^)'^ '^^j^^ . So we come to the following predic- 
tion: decay length for quadrupole particles:Agg ~ ^ — 
is twice smaller than for dipole particles. 




FIG. 2. (Color online) Blue thick line is the border of the 
attraction (inside) and repulsion (outside) zone for parallel 
dipoles in the planar cell from HlOf) . Director no||z. Red thin 

2 

line is the parabola z = Dashed lines make angle tp — 

arccos{-^) with z and are borders of repulsion and attraction 
zone for unlimited nematic. 



B. Interaction in the planar cell with thickness L 

In order to find Green function for this case let's turn 
coordinate system (CS) of the homeotropic cell CS^°™ 
{x,y,z) round the y axis on 7r/2. Then we will have 
(jgpian (j^y^^) with transition matrix yl: x = Ax,x' = 
AS^' so that x = z,y = y, z = —x. Then Ghom{^,^') — 
Ghom{.Ayi,Ayi') = Gp;an(x, x'). Omitting sign ~ we may 
write Green function for planar cell in the C'S'^''^" with 
Ta\\z and x perpendicular to the cell plane [x £ [0, L]): 



. OO OO / 

r^cell I 'i.rn(ip-ip') ■ '^'^^ ■ riTTX 1277 p^ niT Py 

GpianK^^^ ) ^ Y^- 2^ 2^ e >sin—sin^-I^{-^—)Km{^—) (9) 

n— 1 m— — oo 
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where heights x,x' , horizontal projections (0< = \fy^~+^,p> — \f\p~Vz^ , tg^p — ^,tgip' = and p< is less than 
p>. Then taking derivatives brings dipolc-dipolc interaction in the planar cell i/Jrf.p/an = ^^'^^PP'^fj.d'^GpfJ^^: 



U, 



dd^plan 
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{Fi - F2Cos^ip) 



(10) 



where 



Elmr) mrx mrx 
sin sin 
2 L L 

n=l 



Ko{^)+K,C^)\-{nnrcos- 
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„ v-^ , , , mrx mrx , nnp , 



(riTr) sin sin Kn( ) 



When both particle are located in the cen- 
ter of the cell x — x' — ^ we have 

-Er=2,e.eJ"-)'^o(^) and F, 
Y.n=i.odd{™?K2{^). In the limit of smah dis- 
tance p L between particles these functions have 
asymptotics Fi — ^ and F2 — > so that we come 

to the well known resuh Udd = ^^^^^(1 - Scos^ip) 
for p <^ L. In the limit of big distances p > L 
we have |^ = l^:^ + o(-|) with accuracy 5% al- 
ready for p — L. So for p > L it may be written 

Uh,,lan_ = '-^F2{p) • (1 - - SO that 

dipole-dipole interaction is attractive for —(pcl£^<y^c, 

ipc = arccos{^J'\~~^) ~ and is repulsive for 

< < 27r — (if dipoles are parallel each other 
p — p' and vice versa if p = —p' ). In other words 

for p > L dipole-dipole interaction is attractive inside 

2 

parabola z = and is repulsive outside this parabola 



(see Fig. [5]). Practically this parabola z = with 
enough eccuracy confines the attraction and repulsive 
zone even for smaller distances p > 0.3L as it is seen 
from the Fig. [2] All numerical calculations in the paper 
were performed using Mathematica 6, and in all series 
we used summation 

To conclude we have found general approach for de- 
scription of the axial colloidal particles of the size 
0.1pm 10pm in the confined NLC. The decay length 
for dipole interaction is found to be twice more than 
for quadrupole interaction in the homeotropic cell. In 
the planar cell bounding surfaces crucially change attrac- 
tion and repulsion zones for the distances larger than 

2 

Pc = 0.33L where crossover to the parabola z — 
takes place, so that attraction zone is inside this parabola 
and repulsive zone is outside it. This approach has been 
succesfuUy applied as well for the interaction of one par- 
ticle with the one homeotropic and planar wall and for 
interaction between two particles near such wall. That 
results will be published in the upcoming paper [l8| . 
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